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1 Introduction 



In a recent paper we presented numerical results for the nucleon form factors in a fully 
Lorentz covariant model based on the idea that baryons may be viewed as bound states 
of confined constituent diquarks and quarks interacting via quark exchange. The confine- 
ment of the constituents is hereby effectively parametrised. The poles of the corresponding 
propagators are removed by some modifying multiplicative factors which, however, pos- 
sess an essential singularity at time-like infinite momentum. The physical picture behind 
this baryon model is very natural: Diquarks are allowed to "decay" into two quarks, 
one of them recombines with the third quark and forms another diquark. The hope is, 
that for physics relevant at intermediate momentum transfers, most of the complicated 
structure of the baryon may be efficaciously described by assuming strong correlations in 
the quark-quark channel. Thus the notion of diquarks parametrises to some extend un- 
known non-perturbative physics within baryons. In the last years diquarks have not only 
been used in non-perturbative calculations but also in the description of inelastic lepton- 
nucleon scattering, see refs. where various applications of diquarks are extensively 

discussed. Very recently diquark masses have been estimated from lattice measurements 

i- 

The main purpose of our investigations is to formulate a baryon model applicable to the 
intermediate energy region. This is mainly motivated with the advent of a new generation 
of continuous beam facilities like CEBAF at TJNAF, MAMI, EISA, COSY etc. which 
are designed to explore an intermediate region lying between the non-perturbative low- 
energy and the perturbative high-energy regime of QCD. These facilities explore various 
hadron observables to a very high precision. The different existing hadron models are to be 
judged by their ability to predict and to explain these observables in the near future. While 
there exist many models capable to describe pion properties which are strongly dictated 
by chiral symmetry, the detailed structure of all other mesons are still unclear. Although 
quirte a few baryon models, see e.g. 0, |, |^, ||, |, |10|, |lT|, were developed in the last forty 



years, a unified description of baryons within a covariant field theoretical approach, and 
with quarks and gluons as the fundamental degrees of freedom, is still missing. Besides 
being covariant, such a description should include chiral symmetry in its spontaneously 
broken phase and confinement. A covariant model with the correct symmetry pattern 
but without confinement is the Nambu-Jona-Lasinio (NJL) model [|12|, |T3| where quarks 



as fundamental fermion fields interact locally, for recent reviews see [IJ, |16[. There 



baryons appear either as non-topological solitons |T^, ^ |TP[ or as bound states of a 
quark and a diquark ^ |2^, ^ ^ . A hybrid formalism combining the soliton 



with the diquark picture has been developed in ref. |2^, where it turned out, that the 
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soliton background contributes as much to the total binding energy of the nucleon as the 
direct couphng between two quarks and between quark and diquark. In the Global Colour 
Model, a non-local extension of the NJL model, there also exist preliminary studies of 
nucleons as diquark-quark bound states, see e.g. ||2^ and references therein. 



In this paper we extend the investigations reported in ref. There we solved the 
Bethe-Salpeter equation for nucleons in ladder approximation but restricted ourselves to 
scalar diquarks. Employing the solution of the Bethe-Salpeter equation, i.e. the nucleon 
vertex functions, we calculated various form factors like the electromagnetic, the weak and 
the pionic form factor of the nucleon. Despite the fact that this work is a promising start- 
ing point for further investigations, some results, especially for the magnetic moments, 
signalled that the axialvector diquark channel is necessary for a realistic description of 
baryons as bound states of quarks and diquarks. Here we therefore include the axialvec- 
tor diquark channel into the nucleon Bethe-Salpeter equation. We also present results for 
baryons with spin 3/2. 

The paper is organised as follows: In the next section the covariant and confining 
diquark-quark model, as defined in ref. is briefly reviewed. By modifying the prop- 
agators of quarks and diquarks an effective modelling of confinement enters our model. 
In section |^ the Bethe-Salpeter equations for spin-1/2 and spin-3/2 baryons, determin- 
ing their masses and wave functions, are discussed. With an appropriate three spinor 
basis we then construct covariant ansatze for the wave functions which are suitable for 
the numerical solution of the homogeneous integral equations. Section ^ is devoted to 
a discussion of the subtleties concerning the momentum routing in the Bethe-Salpeter 
equation, associated with the ladder approximation. The numerical method is shortly 
described in section ^, and our results for the baryon masses and wave functions are de- 
scribed in section ^. In the last section we finally conclude and give an outlook. Some 
technical issues are deferred to three appendices. 

2 The Covariant and Confining 
Diquark- Quark Model 

Here we briefly recapitulate the definition of the covariant and confining diquark-quark 
model as given in ref. M. Such a model, although without confinement has been also 



considered in ||2^, |30[ . Since the solution of a rigorous relativistic three body equation, a 
Faddeev equation |^T|, is still missing in field theory we follow the path of using diquarks as 
effective degrees of freedom within baryons. They serve as an efficient tool to parametrise 
some of the unknown non-perturbative features of the baryon wave function. As stated 
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in the introduction such approaches have been successful in the past to describe baryons 
in NJL-type models 0, ^ ^ ^ A three-dimensional reduction of the 



fully covariant Bethe-Salpeter model discussed here, can be found in ref. |]32|, where the 
author solved the bound state equations in the Salpeter approximation which ignores all 
retardation effects. 

Having in mind the diquark-quark description arising within the hadronized NJL 
model [0, some basic structures are fixed: To build up a colourless baryon out of a 
diquark and a quark being in the fundamental representation of the colour group SU(3)c, 
diquarks necessarily live in the colour anti-triplet channel. Furthermore, in order to fulfil 
the Pauli principle, scalar diquarks couple via the antisymmetric generators of the flavour 
group, = {p"^^ '^}, axialvector diquarks via the symmetric generators = {p"'^'^"^}, 
respectively ||20| . Our conventions for these flavour matrices are given in appendix 0. 
In the following, however, we deviate from the NJL model as described in Quarks 
and diquarks are treated as elementary but confined particles, see below, whose interac- 
tion, quark exchange, gives rise to quark-diquark correlations strong enough to bind these 
fields to a baryon. 

This can be formalised with the Lagrangian 

C = qA{x){iYd^-i^q)f{-d'^/™^'i)(lA{x) 

+ A^^ix)i-d,d^ - m20/(-9Vm^+)A^(x) 

- \FlMi{-d''l^\,)F^\x) + im^AL,(x)/(-aVm?OA^(x) 

+ e^^^ \g,ql{x)Ci^n%q^{x)A\^^{x) - gl^AaA^YqEixY^s^l^Cqlix)) . (1) 

The quark field is denoted by q{x), the scalar diquark field by A(x) and the axialvector 
diquark field by A'^(x). Their masses are given by the matrices in flavour space mg,mo+ 
and mi+, respectively. For unbroken flavour symmetry they reduce to mg=diag3(mg), 
mo+=diag3(mo+) and mi+=diag6(mi+). In eq. (|ID, capital subscripts denote colour 
quantum numbers. In the kinetic part of the l"*" diquark the non-abelian field strength 
tensor F'"' = d^'A" - ^''A^' + [A'', A^] appears (A^ = A^^tg). The associated self- 
interactions, however, will not be taken into account. The coupling strengths, of Yukawa 
type, between two quarks and the scalar or axialvector diquarks are given by Qs or Qa, 
respectively. In order to take an extended diquark into account, the point-like couplings 
will be supplemented with momentum dependent factors to be defined later. 

An essential ingredient entering our model is the effective parametrisation of confine- 
ment. This is realised by modifying the kinetic terms of the constituents in the Lagrangian 
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(|T]) . Going to Euclidean[| momentum space and using 

= 1 -e-'^(^+"). (2) 

the model quark and diquark propagators, being diagonal in colour and flavour space, are 
given by 

S{p) = %Z^(l-e-'^^'^<^'<), (3) 

D{p) = -^L^(l-e-'^(^'+™o+)/-o+), (4) 
p'^ + mf ■ ^ ' 



0+ 



D 



Due to the numerators, the mass poles of the propagators are effectively screened^]. The 
strength of the screening is described by the parameter d. Whereas the analytic behaviour 
of the quark propagator (0) is similar to the ones obtained from Dyson-Schwinger stud- 
ies of QCD (details about this approach can be found in ref. |^), the justification for 
using confined diquarks is somewhat more involved. As stated above, diquarks are not 
colour singlets and should therefore be confined by the fundamental interaction of QCD. 
Nevertheless most diquark models Q, ^ 0] which describe them as bound states of 



two quarks predict diquarks to be observable particles. Recent studies, however, which 
investigated the system of the quark Dyson-Schwinger equation and the diquark Bethe- 
Salpeter equation beyond the usually employed rainbow-ladder approximation, were able 
to explain why diquarks do not appear in the observable particle spectrum. This mech- 
anism works in the Munczek-Nemirovsky model MT^, as has been shown in ref. 



as 



well as in an extended NJL model PD[. Since both models assume a simplified but quite 



different quark-quark interaction, one may conjecture that the same also holds true imple- 
menting a realistic interaction. For these reasons we use the notion of confined diquarks. 
In our approach we therefore use the propagators defined in eqs. (^, |^) which do not allow 
observable diquarks. In the following it will be seen, that working with confined quarks 
and diquarks leads to the absence of an unphysical quark-diquark threshold for baryons. 

Additionally we introduce the diagonal approximation to the axialvector diquark prop- 
agator by omitting the p^p^ /m\j^ term. 

It is well-known that the full Proca propagator may cause spurious ultraviolet problems. 
These vanish in the special, unitary gauge and have no physical consequences. As such a 

^ We use an Euclidean space formulation with {7^,71^} — 2(5^i,, 7^ = 7^ and pq — Pa'^p- 
^With / = 1 we refer to the propagators as tree level propagators. 




D'-ip) = --r-z;^r' [^\. (6) 
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treatment is, however, beyond the scope of our purely phenomenological investigation we 
have chosen to avoid this problem by using a diagonal approximation to the propagator. 
Its validity is examined more closely in appendix ^. 



3 Bethe-Salpeter Equation for Octet and Decuplet 
Baryons in the Flavour- Symmetric Case 

Using the Lagrangian given in eq. (|T]) we obtain the ladder Bethe-Salpeter equation for 
octet and decuplet baryons. 



X 



^7oip;P) = -2\9a\'s{p,)D'^^{p,)J0^ysi~q)Y^'M-p). (8) 

In this formulation, the Bethe-Salpeter equation involves the matrix-valued Bethe-Salpeter 
wave functions of octet baryons "^siP'jP) and \E'g(p;P) and decuplet baryons \E^^q(p;P), 
respectively, which are projected on positive parity and spin 1/2 or spin 3/2 (see the 
following subsection for their construction). Their flavour part is given by pure octet and 
decuplet states in SU (3) flavour- They depend on the total momentum P of the bound state 
and on the relative momentum p' or p between the two constituents. Mathematically, the 
Bethe-Salpeter equations are equivalent to coupled homogeneous integral equations. The 
numerical method for its solution is presented in sect. ^, see also . 



Although the Lagrangian (|I]) describes a renormalisable diquark-quark theory, at leat 
at one-loop level, the Bethe-Salpeter equations in ladder approximation are formally di- 
vergent in the ultraviolet. The divergence, of course, should then be cured by the wave 
functions. However, to crudely take into account the extended nature of diquarks we 
work, as in ref. |3y, Q, with a finite interaction in momentum space and modify the 



propagator of the exchanged quark according to 

S{q) - S{q) = S{q) . (9) 

This corresponds to a monopole-type form factor. As a consequence, this also removes 
all formal ultraviolet divergencies. Note that we have absorbed the charge conjugation 
matrix C appearing in the Lagrangian (|I|) using the identity C~^S'^{q)C = S{—q). 

Before actually solving the integral equations it is appropriate to find a suitable basis 
for the wave functions appearing in eqs. 
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3.1 Relativistic Three Quark States and their Wave Functions 

When constructing a baryon out of three quarks, its wave function is formally described 
by a spinor of rank three, ipap-y (a, /5, 7 = 1 ... 4). For octet baryons, a convenient basis for 
this multi-spinor can be found by expanding the direct product of the spinors describing 
quarks of flavour h and c with spinor indices (3 and 7 into the complete set of Dirac 
matrices and taking the direct product with a spinor basis of quark a pi|| : 

^Ip^ = {TuU<pCj,)p,. (10) 

C denotes charge conjugation and u = )'(^ )} basis of positive energy 

Dirac spinors in the rest frame describing fermions with spin up (+) and spin down (— ). 
In this representation F and are Dirac matrices to be expanded in the complete set 

{1,75,7^ 757^CT^"}• 

Under a Lorentz transformation S'(A), this wave function transforms according to 

'4^%^ = iSTiA-'P,A-'p)S-')^^,{Su)AS<pS-'C^,)p^. (11) 

The parity transformation is given by 

i^l^^ = 7^(7^F(7^P,p)«)„(7V(7^ Ap)7'C^75)/j„ (12) 

where tt is the intrinsic parity of the baryon and P = {—P,P'^), p = (— p,p^), 7^^ = 
(-7,7'). 

To ensure total spin 1/2 of the wave function, the free Lorentz indices in (^) are to be 
contracted with the independent momenta involved and each covariant will be multiplied 
by a scalar function. 

In the diquark-quark model, only the two choices for (p = {1,7'^75} are taken into 
account which describe the scalar and axialvector diquark. The covariants which lead to 
positive parity of spin- 1/2 octet states, {Tu)a, can be grouped according to scalar and 
axialvector diquark states with the spinor indices P and 7, see table |I], second row. As de- 
scribed above, the free Lorentz indices in these covariants are contracted with the momenta 
P and p and are then to be multiplied with scalar functions. No further symmetrisation 
of the wave function is necessary as this will be provided by the quark exchange. 

Using this classification scheme, the octet wave function can be denoted by 

4 16 



Each of the four covariants F^ which describe the scalar diquark part in the octet baryon 
wave function is multiplied by a scalar function Z^ and likewise there are twelve scalar 
functions Zf multiplying the axialvector diquark covariants Ff'^. 
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Scalar l~)iniiark - ((^^k la 


Axialvpctor l~)iniiaTk - ^k(^^cIo 


Octet 






Decuplet 




(75 rf'' [pTU^])a 


rf e 


{1, ^,^,p^a'^-p,} 


i ^ 


P>'^5PuCT''''Pp, p^^^P.a'^Ppp, l^l^Pu^y'^Ppp] 



Table 1: The Lorentz covariants leading to baryons with positive parity. They are grouped 
according to their diquark content. 



A further reduction of this ansatz by a projection to positive energies is very conve- 



nient. From the expression (0) for the octet wave function we consider only the part with 
spinor index a, multiply it with the adjoint spinor u{P, s) and sum over the spins. This 
leads to a wave function which is, by construction, an eigenfunction of the positive-energy 
projector^] A"*" = (1 + ^). Thus we are led to the following relations between the wave 
functions of eq. ( P^ and the corresponding functions in the Bethe-Salpeter equation (|^: 

The independent covariants Si and Ai are required to be eigenfunctions of A"*" which 
reduces the number of independent scalar functions from sixteen to eight which are now 
denoted by Si, (i=l,2) and Ai, (i=l. . . 6). A convenient representation of these covariants 
suitable for our numerical procedure and for further applications is given by 



S,. 



Si- 


= A+ 


S2-- 


= -pTA^ 


At 


= -^P'^75|^tA+ 


Ai^ 


= P^75A+ 




= PtI^PtA-^ 


A^^ 


= ^P^75A+ 




= 757^A+ - A^, 


A'. 


= ^757^|^tA+-^^ 



(15) 



(16) 



^With a hat we denote normalised four vectors, e.g. PP = 1. In the Euchdean rest frame P — (0, iM) 
this exphcitly reads P — P/iM. Note that aU relative momenta (p, p') are real in Euclidean space as 
they are only needed for spacelike values. 
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where 7^ = 7^^ ~ P^If . Note that the indices have been chosen such that matrices with 
odd indices i = {1, 3, 5} are eigenf unctions to If — iM whereas the ones with even indices 
i = {2, 4, 6} are eigenfunctions to + iM, with eigenvalue in both cases. 

In the rest frame of the bound state, P = {0,iM), the ansatz for the matrix valued 
nucleon wave function which we used for further numerical processing reads 



V *8(p;P) ) 



ISi 

ii<^p)S2 

'^i^P)Ai 
IA2 
ip{ap)A3 + {a- X p){(Tp)A5 
ipA, + l{(Txp)A, 



(17) 



The scalar quantities Si and Ai depend on p = y/p^p^ and z = cos ip = P ■ p for each 
value of the bound state mass. As expected from the properties of Si and Ai (see eqs. 
(15) and (|l^)), upper components have odd indices, lower components have even indices. 

The strategy shown above for the octet baryon can also be applied to decuplet baryons 
which have spin 3/2. Projection of the tri-spinor wave function ipa/S'y onto total spin 3/2 
may be achieved by expanding the piece associated with quark a in terms of Rarita- 
Schwinger spinors u^^: 

= {Tu,U<l>C)p,, (18) 

r is chosen such that all Lorentz indices are contracted with the momenta P and p. 
Note that the Rarita-Schwinger constraints demand u^'j^ = u^P^ = and that in the 
diquark-quark model the diquark part of decuplet states is made of axialvecor diquarks 
only, (p = C'j'^. 

The covariants F which are left for decuplet states after applying the above restrictions 
and the positive parity constraint can be found in table |l|, third row. The covariants for 
octet and decuplet states are closely related. Due to the Rarita-Schwinger constraints the 
Lorentz index of must be contracted either with the Lorentz index of the axialvector 
diquark (7^6*) or with the transversal relative momentum pj.=p^ — P^{p ■ P). For the 
first choice, we have the same four covariants as for the octet state with scalar diquark 
correlations. For the second choice, we obtain the twelve covariants as for the octet state 
with axialvector diquark correlations, however, multiplied by 75 to ensure positive parity. 
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The decuplet wave functions reads now: 

= (Tr^'o«'')a(7.C),, (19) 

4 16 
1=1 1=5 

where the denote the sixteen scalar functions for the decuplet state. 

Taking from eq. ([19|) only the part with spinor index a and projecting it to positive 
energies leads to 

+3/2 

Y: T^o (p; s)u\P, s) = T^^{p- P)^^' = vi/^o^(p; P) = ^^,^{p; P)P'\ (20) 

s=-3/2 

Therefore the projected wave function (which besides being a 4x4-matrix has tensor 
character) is determined by the condition 

vI/^o'(p;P) = <(p;P)P^ (21) 

which requires \E'io to be an eigenfunction of the Rarita-Schwinger projector. Here, the 
explicit expression of the Euclidean Rarita-Schwinger projector is given by 



p/^^ := A+ ( -^t^^ + -^t'Y TT- + ] =■■ A+A''^ (22) 

^ 3 ' ' 3 M2 3 M ' ^ ' 



The most general form which fulfills condition (^Tj) is 

^^o'(p; P) = AcS.A+A^'^ + zE,75^rA+p^A^'^ (23) 

which requires the covariants again to be eigenfunctions of A"*", as a consequence we 
may use the same choice for the and Si as in the octet case and are left with eight 
independent scalar functions. 

In the rest frame of the bound state, the decuplet wave functions are then denoted by 



3 

' 3' 

■ ^{pi -l{ap)a^)E2 0\ 
i{cTp){p ~\{cTp)a^)E, o) 

and all other components of "^iq {p;P) vanish. The appearance of 75 in eq. ( P^D has 
interchanged upper and lower components as compared to the octet case. 
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Instead of working with the Bethe-Salpeter wave functions one may alternatively use 
the Bethe-Salpeter vertex functions obtained by amputating the external quark and di- 
quark propagators from the wave function: 

Mp;P) = s-\pa)D-\p,)^s{p;P), (25) 

mp;P) = S-\pa){D-'r{p,)%{p;P), (26) 
$^o'(p;P) = S~\p,){D-'r{PbmiP;P), (27) 

Substituting the wave functions by the vertex functions in the Bethe-Salpeter equations 
(|^^ leads to a reformulation of the bound state equations which is sometimes more 
convenient. For example in ref. |3y, |1| the equations containing only scalar diquarks have 



been solved in a form which includes the vertex function explicitly. 

3.2 Orbital Angular Momentum and Spin of the Bethe-Salpeter 
Wave Functions in the Rest Frame of the Bound State 

Whereas the choice of the covariants in (1^.16) which build the octet and decuplet wave 
functions is well suited for numerical computation and further covariant calculations, their 
physical interpretation is not obvious. 

In general, covariant wave functions possess only the mass of the bound state M and 
its total angular momentum J as good quantum numbers. In the rest frame of the bound 
state, however, the wave functions can be written as a sum of tri-spinors each possessing 
definite orbital angular momentum and spin, thus allowing a direct interpretation of the 
different components. These tri-spinors are linear combinations of the covariants Si and 
Ai which have been constructed in the previous subsection, multiplied by the respective 
Dirac matrices (75C)/37 and {'y^C)i3^ denoting the diquark content, respectively. 

In the rest frame the Pauli-Lubanski operator for a tri-spinor is given by 

W' = \e^,kiy\ (28) 
whose square characterises the total angular momentum 

iyWVa/37 = ^(^+l)^a/37- (29) 

Here, is the tri-spinor wave function with positive parity and positive energy. The 
tensor is the sum of an orbital part, ^ and a spin part, S*-^^, which read 

2(5"-' )aa',/3/3',77' = io^ )aa' ® ^/3/3' ® (^77' + ^aa' ® )fifi' ® (^77' + 

® 5/3/3' O (0^-'''=)77'' (31) 
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such that C^'^ = + ^S^^. Obviously, is proportional to the unit matrix in Dirac 
space. The definition of a^'^ := — 1[7'^,7'^] differs by a minus sign from its Minkowski 
counterpart. The tensors L and S are written as a sum over the respective tensors for 
each of the three constituent quarks which are labelled a = 1 ... 3 and with respective 
Dirac indices aa', 77'. 

With the definition of the spin matrix = \eijk(J^^ the Pauli-Lubanski operator reads 



■ d 

V = {-i)eijkP'- 



1 



i 



77 



13/3' ^ fJ^y 



Syy'-\- 



'77 y ' 



(32) 
(33) 

(34) 



where we have already introduced the relative momentum p between quark and diquark 
via a canonical transformation: 



P = p^+p^+p^, p = ri{p^ + p^) - {I - r])p^ , P=-^{P^-P^)- 



(35) 



Assuming a pointlike diquark, the relative momentum between quark 1 and 2, p' , vanishes 
and the only contribution to the orbital angular momentum stems from p. now takes 
the form 



2(L ■ S) 



Jao',/3/3',77' 



d 

-ejjfcF'^ ((S*)a„/ (g) 6(313' ® (^77' + 

® [(S%;3' ® ^77' + 6,3(3' ® (S%y 
^ (9 6aa' ® 5/3/3' ® ^77' + 2 5aa' (g) (S*)/3/3' ® ^77' + 5/3/3' ® (S*)77' 
+2 5„«' ® (S')/3/3' ® (S%y 



(36) 
(37) 

(38) 

+ 
(39) 



When applying to the wave functions we first note that a scalar function does not 
contribute to the angular momentum, e.g. 



Si{p\ iMpi) = 2L-S Si{p\ iMpi) = 0. 



(40) 



In table ^ we express all terms of the octet wave function in terms of eigenstates of 
and S^. To this end we define the following linear combinations of matrices: 



3) 



12 



ipaB'y ^^^^ frame 



52m(75C) 



X 



(75C) 













'-(f -\a\cTp))x. 



eigenvalue 
l{l + 1) of 




2 

2 



2 

6 
2 



eigenvalue 

s{s + l)of 



15 
4 

15 
4 



Table 2: Classification of the components of the octet wave function in terms of eigenstates 
of and in the rest frame of the bound state. 



= —A^e + lAl (41) 

The eigenvalue of is | for all terms, of course. In the table, x = {x^? X~} denotes an 
arbitrary Pauli two-component spinor which is the positive energy basis for quark a with 
Dirac index a. To derive this, the following relations between Dirac matrices have proven 
to be useful, 

S^'(75C) + (75C)(S^7 = (42) 
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ip^a^ in the rest frame 



eigenvalue 

/(/ + 1) of L2 



eigenvalue 

s(s + l)of 



V(o-p)(px)/ 



'(y\crp){px) 







(rc) 



iYC) 




2 

2 
6 

6 
2 

6 
12 



15 
4 

15 
4 



15 
4 

15 
4 



Table 3: Classification of the components of the decuplet wave function in terms of eigenstates 
of and in the rest frame of the bound state. 



E^(75C)(E^7 = -3(75^), 



(44) 



Three covariants can be regarded as "s-wave" components, <Si, A2 and Bi, and we 
expect the corresponding scalar functions to dominate the wave function decomposition. 
The other "p-, d-wave" components represent all remaining possibilities of combining 
orbital angular momentum between quark and diquark and the joint spin of axialvector 
diquark and quark to total spin 1/2. In this sense the description is closed. 

The individual terms of the decuplet wave function can be classified accordingly with 
the help of the linear combinations: 



Si, 



2i 
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Sr = {-jl,Af^ + ll,A^ P^T - 1^,6^" ■ (46) 

The result can be found in table |. Due to the spin projection by use of the Rarita- 

Schwinger projector the eigenvalue of is ^. = {x^\X~^} (^=1)2,3) denotes a 

two- component vector-spinor which survives the spin-3/2 projection in the rest frame. 

The Rarita-Schwinger constraints reduce to crx = 0. 

For the decuplet only one s-wave component exists, Pi. Again these covariants exhaust 

all possible couplings of spin and orbital angular momentum, note that even an orbital 

angular momentum 1=3 contributes to a spin-3/2 state, Sq. A contribution of an 1=1 

/ \ 

state with scalar diquark, such as (75^"), is forbidden by the Pauli principle for 

V (PX) J 

pure decuplet states, but may admix in the case of broken flavour symmetry. 



4 Ladder approximation to the Bethe-Salpeter equa- 
tion 

In our preceding publication we used a momentum routing in the Bethe-Salpeter equa- 
tion, where the interaction (propagator of the exchanged quark) is manifestly independent 
of the total momentum P of the baryon, i.e. q = —p — p'. A P-independent kernel is 
desirable for the following reasons: As will be discussed in the next section, it reduces 
the numerical work when solving the integral equation quite drastically. Furthermore, 



the canonical normalisation condition for the Bethe-Salpeter wave function |^ becomes 
much more involved when the interaction is P-dependent. Therefore, when solving a 
Bethe-Salpeter equation in ladder approximation one usually tries to find a momentum 
routing having this very convenient feature. In our equations, however, where quark and 
diquark interact through quark exchange, which changes the identity of the particles af- 
ter each interaction (a diquark becomes a quark and a quark becomes a diquark after 
the quark exchange), demanding a P-independent interaction defines a "modified ladder 
approximation" with results deviating from other momentum routings. Note, that this 
problem does not occur in the Bethe-Salpeter approach if the interaction, e.g. meson ex- 
change in the two nucleon system or gluon exchange between quark and antiquark, does 
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Pa=^P-P 



Pa=^P-P 





Pt,=(l-Tl)P+p 



pJ=(l-Tl)P-p' Py=(l-Tl)P+p 



Figure 1: Momentum routing defining the modified ladder approximation. 

not change the identity of the constituents. Almost all Bethe-Salpeter equations treated in 
the literature are of this kind. For a comprehensive review of the existing Bethe-Salpeter 



literature see In order to clarify the ambiguity arising in the diquark-quark Bethe- 
Salpeter equation, which came up during our investigations, we discuss here two possible 
momentum routings and in section 6 and appendix P the corresponding variations of the 
results. 

For sake of clarity of this discussion we temporarily restrict ourselves to the part of the 
octet equation which involves only scalar diquarks, i.e. eq. (|^ with ga = 0. Furthermore, 
it is more transparent to work with the Bethe-Salpeter equation which involves the vertex 
function (see eq. (P^D): 



\9s 



(27r) 



:l,S{-qh,S{p'jD{pi)^s{p';P). 



(47) 



Repeating the procedure discussed in section ^Tl] for the vertex instead of the wave function 
leads to an ansatz [BOl PI 



Mp; p) 



lSi{p,z) 

MCTP)S2{P,Z) 



(45 



Note that in the lower component of this spinor-like object the spin is orientated along 
the spatial part of the relative momentum. In the Bethe-Salpeter equation (|47|) the 
momentum of the constituent quark is given by p'a = p' + vP^ the momentum of the 
scalar di quark by p'^ = —p' + {1 — ri)P. P denotes the total momentum of the nucleon 
and p the relative momentum between quark and diquark. The Mandelstam parameter 
7] describes how P is partitioned to quark and diquark. The fact that the eigenvalues 
of the Bethe-Salpeter equation are independent of r/ is a direct consequence of Lorentz 



covariance. 



p', the interaction 



If the momentum of the exchanged quark is chosen as g = —p 
kernel is by construction independent of the total momentum. Since the quark exchange 
transforms a quark to a diquark and vice versa, after the interaction, the relative momen- 
tum between these particles becomes (1 — ri)pa — rjph = —p, see fig. |T] for this momentum 
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routing. We now demand ^s{~P', P) = ^sip, ~z,P), i.e. the orientation of the nucleon 
spin should not depend on the orientation of p as required by the Dirac decomposition 
before projection onto positive energies. Then the vertex function appearing on the left 
hand side in eq. (^7]) reads 



lSi{p,-z) 
Map)S2{p,-z) 



(49) 



The definition of the momentum q together with eq. (|49|) is what we call "modified ladder 
approximation" which leads to satisfying results for various nucleon form factors in the 
weak binding regime, as can be seen in ref. Using such a prescription we have been 



able to reproduce the results of ref. ||30 



Another possible momentum routing allows P to flow through the quark exchange, 
q = —p — p' + (1 — '2ri)P, see fig. 0. When taking not into account, that quark and diquark 
change their role after exchanging a quark, that means attributing the index "a" to the 
quark and index "b" to the diquark, the relative momentum after the interaction is given 
by p. Correspondingly, the vertex function on the left hand side in eq. (^7]) is denoted by 



$8(p; P) 



lS,{p,z) 
^iap') S2{p, z) 



(50) 
(51) 



We name this "direct ladder approach". Although one might expect that both choices 
of the momentum routing lead to the same physical results, we found that this is actu- 
ally not the case. In the next sections and appendix P we report on calculations using 
both momentum routings, see also for further results obtained in the direct ladder 
approach. We found that physical results slightly differ for these two choices. We will 
show, however, that both possibilities are manifestly Lorentz covariant, i.e. the eigen- 
values do not depend on the Mandelstam parameter rj. Furthermore, for weak binding 
the eigenvalues almost coincide. Given that the ladder approximation is only reliable for 
weak binding we conclude that both methods are of similar validity. 



Pj=(l-Tl)P-p 




P(,Kl-Tl)P-p 




p=riP+P 



q=-p-p'+(l-2Tl)P 



pi=(l-Tl)P-p' Pa=nP+P 



Figure 2: Momentum routing defining the direct ladder approximation. 
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Although during this discussion we restricted ourselves for clarity to scalar diquarks 
and to the Bethe-Salpeter equation involving the vertex function, it can be extended in a 
straightforward way to the complete equation including the axialvector diquark channel, 
the decuplet equation and also to the equations (^,||) involving the wave functions in an 
obvious way. The reported numerical results are always given for the full problem. 



5 Numerical method 



For the numerical solutions of the Bethe-Salpeter equations we developed an iterative 
hybrid algorithm, which allows a very efficient and fast computation. A description of 
our numerical method can be found in a forthcoming publication which presents this 



algorithm in all details Thus we will focus here on the main steps only. 

We solve the Bethe-Salpeter equation as a system of equations for the wave function \1/ 
and the vertex function $, see eq. Both wave and vertex function can be expanded 

in the rest frame according to eq. (^, Although there are 10 equations for the eight 
octet functions, stemming from 2 equations in the scalar diquark channel and 2x4 equa- 
tions in the axialvector diquark channel, we confirmed that two of them are redundant. 
The decuplet system yields 2x12 equations for eight scalar functions due to the tensor 
character of the wave function which reduce again to eight independent equations. This 
especially underlines the necessity of including the sub dominant amplitudes describing 
orbital angular momentum to keep the system closed. 

We expand the scalar functions (amplitudes) S, A and D, E defined in (0,^) into 
Chebyshev polynomials of the second kind. 



n=0 

oo 



(52) 
(53) 



n=0 



where amplitudes with a hat, i.e. Yi, belong to the vertex function and the ones with- 
out a hat to the wave function. Here we use the generic label Yi according to Yi = 

{Si, 5*2, . . . , ^ejoctei) {-Dl, D2, . . . , EQ}(iecuplet- 

Throughout the calculation we work with usual hyper spherical coordinates. 



fp'A 

P's 

\p'J 



p 



[ sin ip' sin 9' sin 0' ^ 
sin sin 9' cos 0' 
sin ip' cos 9' 

\ cos Tp' J 



cosip' = P ■ p' . 



(54) 
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We are free to choose the spatial part of the relative momentum p appearing on the 
l.h.s. of the Bethe-Salpeter equation (|^. Without loss of generality, we select = 

{0,0, pVT^,pz). 

We expand quark and diquark propagators into Chebyshev polynomials as well and 
project the Bethe-Salpeter equation onto the Chebyshev moments of the amplitudes, 
Yi"^, Y^. Note that in the chosen Lorentz frame with the relative momentum p parallel 
to the third axis this is especially easy because in this case the amplitudes as given in 
(0) do not mix. The integration necessary to generate the kernel will be performed in 
hyperspherical coordinates. The integrations over 0' and 9' are done analytically, and the 
remaining two over z' and z (due to the projection) numerically. 

The final equation suitable for iteration or diagonalisation reads: 

j=l n=0 

ypiPi^) = E E / dp'p''H^^{pi„pl)Y,"^{p'J. (56) 

k=l m=0 

Here, P/^" is the propagator matrix with amplitude indices i,j and Chebyshev mo- 
ment indices m,n and Hj'^{p,p') is the matrix of the quark exchange kernel given on a 
momentum grid {pii,p'i,J with respective amplitude indices j,k and Chebyshev moment 
indices n, m. The sum over amplitude indices runs from 1. . . 8 for both octet and decuplet. 
rrimax and rimax denote the highest Chebyshev polynomial considered in the expansion of 
the vertex and the wave function amplitudes, respectively. The kernel includes also the 
flavour factors and the ratio Qa/ Qs (for spin-1/2 baryons). It is only in the modified lad- 
der approximation that the integral kernel does not include an explicit dependence on the 
bound state mass M for all values of rj, which makes a fast determination of M for given 
couplings feasible. We refer to g as the eigenvalue of the Bethe-Salpeter equation which 
is Qs for spin-1/2 and ga for spin-3/2 baryons. 

In the following we will investigate the convergence properties of the expansion in 
Chebyshev polynomials. Furthermore, the independence of the eigenvalue from the Man- 
delstam parameter r] will be discussed. The chosen example is the octet equation (|^) with 
the axialvector diquark propagator diagonal in Lorentz indices and a diquark size factor 
of the monopole type. 

For all calculations reported, a momentum mesh size of 20x20 and the inclusion 
of only zeroth, first and second Chebyshev moment for both vertex and wave function 
amplitudes is sufficient for determining the eigenvalue up to 10~^ precision (see table 
This is also reflected in the magnitudes of the vertex amplitudes Yi. Going up one 
Chebyshev moment suppresses the amplitudes by almost one order of magnitude. The 
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Confining Propagators, d=10 




^max 





1 


2 


3 



















9.4362 


9.1533 


9.1417 


9.1420 


1 




9.7480 


9.2277 


9.2260 


9.2265 


2 




9.7568 


9.2276 


9.1992 


9.1988 


3 




9.7568 


9.2277 


9.1994 


9.1990 


Confining Propagators, d=l 




^max 





1 


2 


3 



















11.3797 


10.9507 


10.9544 


10.9547 


1 




11.8555 


11.1401 


11.1704 


11.1711 


2 




11.8656 


11.1503 


11.1443 


11.1445 


3 




11.8656 


11.1504 


11.1444 


11.1446 



Table 4: Octet Bethe-Salpeter Equation - eigenvalues Qs for various rrimax and rimax, the 
maximum Chebyshev momenta of vertex and wave function. Parameters: mg=mQ+=mi+ , 
M=1.5mg, A=2mg, ga/gs=0-^, i]=0.5, momentum grid size np=20. 

wave function amplitudes Yi converge somewhat slower, see also sect. |6.3| . For a large 
pole screening factor, d > 5, and a weak binding situation the even Chebyshev momenta 
are more pronounced than the odd ones. 

Increasing the couphng constant Qa in the nucleon equation always lowers the eigen- 
value, hence the axialvector diquark enhances binding. 

For tree level propagators the choice of the Mandelstam parameter rj is limited to the 
values rj & [1 — min{mo+, mi+}/M, mg/M] to avoid singularities in the propagators. This 
restriction should not apply for confining propagators. To demonstrate this, we choose 
the extreme case mo+=mi+=2mq and the tree level threshold M=3mq0. Then, tree level 
propagators limit the choice to ?7=l/3, as opposed to confining propagators. For a small 
pole screening factor d=l we could vary r/ in a wide range without affecting the eigenvalue, 
for d=10 due to numerical instabilities the invariance region depends on the constituent 
masses, displayed in figure ^ in the left panel. It can be seen that the poles in the tree level 
propagators are effectively screened although the invariance region for d=10 is restricted 
to 0.2. ..0.42. 

Whereas the latter results were obtained in the modified ladder approximation, the 

*For similar choices of quark and diquark masses we calculate the octet and decuplet masses, see sect. 
|. The masses of the decuplet baryons are close to the tree level threshold. 
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Figure 3: Eigenvalue vs. Mandelstam parameter rj, left: two pole screening factors compared, 
right: modified and direct ladder approach compared. Dotted lines denote the location of the 
propagator poles. 

check of ?7-independence for the direct ladder approach requires the total momentum P to 
appear in the kernel: q = —p—p' + {l — 2ri)P, which slows down the numerics considerably. 
In the case of tree level propagators, the right panel of figure |^ shows the invariance of 
the eigenvalue for both momentum routings. 

The crucial advantage of using confining propagators can be seen even more clearly 
in figure ^ Here the function gs{M) decreases rapidly for tree level propagators near 
threshold while the corresponding function for confining propagators runs smoothly over 
the "pseudo" -threshold. We furthermore observe that the even Chebyshev momenta of 
the Bethe-Salpeter wave functions for tree level propagators close to threshold become 
squeezed in the low momentum domain while the corresponding odd momenta are sup- 
pressed. No such effect is present for the Bethe-Salpeter wave functions obtained with 
confining propagators. 

All the numerical features described also apply to the decuplet equation in the weak 
binding regime (M > + mo+,i+))- For lower bound state masses the convergence 

becomes worse and for M < + mo+,i+) our method failed in finding a real eigenvalue 
for the coupling Qa- In this mass region an effect well known from the Cutkosky model 
is visible: two states described by their functions ga{M) collide and form a complex 
conjugate pair of eigenvalues This shortcoming of the ladder approximation 

provides us with an upper limit of the coupling constant Qa for which the approximation 
is valid. However, the ga needed to describe the baryon decuplet is well below this critical 
value. 
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Eigenvalue vs. Bound State Mass 


m,=ni„*=m/, A=2m^, gjg 


,=0.5 




- 
- 


Tree level prop. 


\ 

\ 


Confining prop., d=1 


X 

X 


Confining prop., d=10 


X 







1.5 1.6 1.7 1.8 1.9 2.0 2.1 2.2 2.3 



Figure 4: Eigenvalues vs. bound state mass M. Absence of threshold effects for confining 
propagators. 

6 Results: Masses and Selected Wave Functions 



6.1 Octet and Decuplet Masses 

In our approach the strange quark constituent mass rris is the only source of flavour 
symmetry breaking. Isospin is assumed to be conserved. The equations describing octet 
and decuplet baryons have been derived under the premises of flavour and spin conser- 
vation, i.e. only wave function components with same spin and flavour content couple. 
The flavour structure of the eight equations describing A^, A, S, S, A, S*, S* and Q can be 
found in appendix 0. 

In order to limit the number of parameters we assume the scalar and axialvector 



diquark masses to be equal. Furthermore, we choose them to be mg+ ^+ = ^{nif + 
where fg G {uu, us, ss} is the flavour content of the diquark. We denote ^ as diquark 
mass parameter. Assuming ^ G [0, 1] is obviously natural. 

With these choices the model has the following parameters: two constituent quark 
masses m„ = md,ms, the pole screening factor d, the diquark size factor A, the diquark 
mass parameter ^ and the couplings ga,gs- 

We do not try to make a thorough fit onto the baryon masses in our parameter space. 
Inspired by our results for the form factors of proton and neutron we assume first 
d=10. This results in only slight modifications of the propagators compared to the tree 
level ones for spacelike momenta. Furthermore, we choose A=l GeV ~ 2mu and ^=1 to 
stay close to calculations of octet and decuplet masses in NJL-diquark-models as done 
in and These authors, however, used just a static approximation to the Bethe- 
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exp. 


mod. 


ladder 


dir. ladder 






I 


II 


III 


d 




10 


1 


1 


A (GeV) 




1 


1 


1 


mu (GeV) 




0.5 


0.5 


0.5 


rus (GeV) 




0.65 


0.63 


0.63 






1 


0.73 


0.73 


9a 




10.35 


10.92 


10.05 


9s 




9.43 


8.06 


7.34 


Ma (GeV) 


1.116 


1.123 


1.130 


1.133 


Ms (GeV) 


1.193 


1.134 


1.137 


1.140 


Mh (GeV) 


1.315 


1.307 


1.319 


1.319 


Ms* (GeV) 


1.384 


1.373 


1.372 


1.380 


Ms* (GeV) 


1.530 


1.545 


1.548 


1.516 


Mn (GeV) 


1.672 


1.692 


1.697 


1.665 






0.0028 


0.0028 


0.0021 



Table 5: Octet and decuplet masses obtained with the maximum order in Chebyshev polyno- 
mials 'rnmax=nmax=^ and momentum grid size np=20. 

Salpeter equation. Note that this part of our calculation has been done only in modified 
ladder approximation since we were forced to choose t] <0.4 (cf. figure |]), and solely for 
ri=0.5 the integral kernel does not depend on the baryon mass in the direct approach, 
which allows a fast determination of M. We then fix the couplings Qa and Qs by the 
nucleon and delta mass. The ratio ga/9s is quite independent from the ratio M^/rUg and 
varies weakly with d. Finally we vary just the two constituent quark masses to obtain the 
other six hyperon masses reasonably close to their experimental values. 

Secondly, we explore the case d=l and now additionally allow the diquark mass factor 
^ to vary. In |^ we showed that this strong screening of the propagator poles leads 
to overestimated nucleon e.m. form factors for high Q^. In this case we also compare 
results for the two momentum routings, modified and direct ladder approach. In both 
calculations r] is set to 0.5. 

The results are given in table |^. The quality of the results may be read from the 
quantity = J2^=i ^^^''''^ . Note that the mean deviation of the calculated masses 

from the experimental ones is of the order of half a per cent or less. Column I shows 
the results for d=10. A remarkable feature is the large constituent quark mass 171^=500 
MeV, whereas the constituent mass difference between strange and up quark is 150 MeV, 



23 



a commonly used value. In the following two subsections we will show for this set of 
parameters the wave function amplitudes for the nucleon and the delta and discuss the 
contribution of the various components with respect to their orbital angular momentum 
and exemplify flavour symmetry breaking effects on the wave and vertex function of the 
A-hyperon. 

The next columns, columns II and III, show the results for d=l with the other param- 
eters chosen to give baryon masses to the same level of accuracy as before. In this case 
quark and diquark propagators are strongly modified for spacelike momenta. Quark and 
diquark "masses" begin to loose their meaning which we usually attach to them. This is 
reflected in our solutions in a rather small diquark mass parameter which renders con- 
stituent quarks and diquarks roughly equal in their masses. Direct and modified ladder 
approach give approximately the same results in the latter case. 

In all cases the Gell-Mann-Okubo mass formula for the octet is fulfilled with an inac- 
curacy of less than 0.5%. It should be emphasized that the mass sphtting between octet 
and decuplet is exclusively provided by the quark exchange with the couphng strengths gg 
and ga and is not due to a heavier axialvector diquark as assumed in other diquark-quark 



models, e.g. in ET 



6.2 Vertex Functions for the A-hyperon 

The A-hyperon appears to be of special interest. First, its measured polarisation asym- 
metry in the process p~f —>■ K^A could provide a stringent test for the diquark-quark 



model. As discussed in [^, there are only scalar diquarks involved in this process and in 
the following we will concentrate on the scalar diquark part of the vertex function. 

Secondly, broken S'f/(3)-fiavour symmetry induces a component of the total antisym- 
metric flavour singlet [{su)d + {ud)s + ((is)^]^ into wave and vertex function (see ap- 
pendix As the flavour singlet is only composed of scalar diquarks and quarks, this 
generates two additional scalar amplitudes Si^singiet and S2^singiet besides the usual two 
from the octet A state {su)d — \/2{ud)s + {ds)u . The vector part of the vertex func- 
tion remains unchanged in flavour space, ^ {[suYd — [dsYu). 

These scalar amplitudes for the vertex function are depicted in figure ^ for the pa- 
rameter set of column I in table |[ The S'l^smg^ercomponent is suppressed against Si^octet 
by two orders of magnitude. However, the purely relativistic S2^singiet- component is only 
5 times smaller than its octet counterpart. Observing that the S'l-component usually 



^With round brackets we denote scalar diquarks, e.g. [us) — us — su. Axialvector diquarks we denote 
as [its] = us + su. 

^In nonrelativistic quark models with SU{&) symmetry such a component is forbidden by the Pauli 
principle. 
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Figure 5: Scalar amplitudes of the vertex function for the A-hyperon, normalised to 
S'°orfet(Pi)=l with Pi being the first point on the momentum mesh. The used parameters are 
m„=0.5 GeV, m,=0.65 GeV, ^=1, r/=0.33, rf=10 and A=l GeV. 

contributes the major part to observables as demonstrated in we can safely regard the 
A-hyperon as an almost pure octet state in flavour space. 
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6.3 Wave Functions for Nucleon and Delta 



In this subsection we present Bethe-Salpeter wave functions for nucleon and delta, using 
the parameter set of column I in table |^. They are normalised to Si{pi)=l (nucleon) or 
Dl{pi)=l (delta) where pi is the smallest point of the momentum mesh. The amplitudes 
represent the strengths of the (L^, S^) eigenf unctions given in sec. |3.2| and are simple 
linear combinations out of the amplitudes defined in (0,^). 

As already mentioned, the convergence of the wave function amplitudes in terms of 
Chebyshev polynomials is somewhat slower than for the vertex amplitudes: Second and 
zeroth Chebyshev moment of the amplitudes differ by less than one order of magnitude. 
All wave function amplitudes are concentrated to four-momenta p < 0.6 GeV. 

In figure |^, the nucleon amplitudes with even orbital angular momentum / appear 
in the left row: These are the three s-waves describing (i) scalar diquark and quark, 
(ii,iii) axialvector diquark and quark oppositely aligned to give spin 1/2. There are two 
axialvector diquark components due to the virtual time component of the latter. The 
scalar diquark component is the most important but the other s-waves enhance binding 
by approximately 30 %. The fourth "non-relativistic" component is a strongly suppressed 
(i-wave with quark and axialvector diquark aligned to give spin 3/2. The lower components 
depicted on the right side can be understood as the admixture of negative-energy spinors 
to the proton wave function and contribute approximately 10 % to the binding energy. 

The delta amplitudes in fig.0 have also been arranged into "non-relativistic" (left row) 
and "relativistic" components (right row). As expected, the only s-wave dominates the 
decomposition, but the relativistic p-wave components, which act repulsively, increase the 
eigenvalue by approximately 20 % and are thus non-negligible. Again the component with 
the highest orbital angular momentum, 1=3, is highly suppressed. 

No clear effects of flavour symmetry breaking, as in terms of the width in momentum 
space, can be detected for the other octet wave functions when compared to the nucleon. 
For the vertex functions (where the zeroth Chebyshev moment dominates), we observed 
a limited increase of the width in momentum space with increasing number of s-quarks 
(10% difference between and S). The width is hereby defined as the absolute value of 
the relative momentum p where the zeroth Chebyshev moment of an s-wave amplitude 
reaches half the maximum value. 
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Figure 6: Scalar and axial vector amplitudes of the nucleon wave function with parameters 
given by m„=0.5 GeV, ^=1, r/=0.33, d^lO and A=l GeV. 
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Figure 7: Axial vector amplitudes of the delta wave function with parameters given by m. 
GeV, e=l, d=10 and A=l GeV. 
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7 Conclusions and Outlook 



In this paper we extended the covariant and confining diquark-quark model of ref. [|^ 
by including axialvector diquarks in the description of baryons within the Bethe-Salpeter 
approach. Thus we were able to calculate the octet baryons masses as well as the decuplet 
masses and wave functions. We implemented confinement via an effective parametrisation 
of the constituent propagators and demonstrated the existence of bound states beyond 
the pseudo-threshold. 

We decomposed octet and decuplet vertex and wave functions in the Dirac and Lorentz 
algebra, obtaining 8 scalar functions, respectively, which we computed numerically. Two 
approximations to the Bethe-Salpeter equation were discussed and compared: The direct 
ladder approximation, and a modified ladder approximation with a special momentum 
routing leading to a kernel independent of the bound state mass. 

In order to fix the parameters of the model preliminarily, we computed the masses 
of octet and decuplet for broken SU{3) flavour with isopin conserved. With the scalar 
and axialvector diquark mass assumed to be equal, the octet-decuplet mass splitting is 
a result of the different effects of the coupling constants in the scalar and axialvector 
diquark channel. For a parameter set which fits the octet and decuplet masses well, we 
computed vertex amplitudes and wave functions for all octet and decuplet baryons. The 
wave functions for baryons with different strangeness content but same spin differ mostly 
due to the different flavour Clebsch-Gordan coefficients, the respective scalar functions 
being very similar. Therefore we presented only the wave functions for the nucleon and 
delta out of this data. The decomposition of the wave functions in the rest frame of 
the bound state in terms of spin and orbital angular momentum eigenstates revealed an 
s-wave dominance in all ground state baryons stemming from both scalar and axialvector 
diquark contributions. The p-wave contributions sneaking in via the lower components of 
the spinors are of greater importance for decuplet baryons than for octet baryons. 

Due to its special role among the other baryons, we investigated the A-hyperon in more 
detail and discussed its vertex amplitudes. In our approach, the A acquires a small flavour 
singlet admixture which is absent in SU (6) symmetric nonrelativistic quark models. 

This work, together with our preceding paper ||l|], provides a sound basis for further 
applications of this approach to baryon phenomenology. The calculated amplitudes which 
encode the non-trivial information of a baryon as a diquark-quark bound state serve 
as a necessary input for the calculation of various observables. A calculation of the 
electromagnetic form factors of the octet and decuplet baryons is hereby the next task. 
Such an investigation will furthermore help to flx some of the parameters. Our aim is 
however to apply the covariant and conflning diquark-quark model to processes, which 
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are far less understood. Especially, the reactions p7 KA and pp pAK, currently 
measured at ELSA and COSY, respectively, will serve as a stringent test of our approach. 
Additionally, we plan to get further insight into this picture of baryons by computing 



structure functions for a spectator model which includes the axialvector diquark. 
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A Octet and Decuplet Equations 

The symmetric and antisymmetric flavour matrices can be written down as: 



,a 
I A 




(57) 



(58) 



By these conventions flavour antisymmetric diquarks are (ud) := ud — du, etc. and 
flavour symmetric diquarks \piuu, [ud] := ud ^ du, etc. The flavour wave functions of 
octet and decuplet states do not decouple once the s-quark breaks the symmetry. The 
Bethe-Salpeter equation (^ still describes nucleons (isospin is assumed to be conserved), 
and equation (^ still refers to A and Q which possess only single-component flavour wave 
functions [uu^^u resp. [ss]'^''s. 

We use the following abbreviations to give short-hand Bethe-Salpeter equations for 
the remaining baryons: 



30 



{ab)c resp. [a6]^c 


: octet flavour wave functions witli diquark flavour 
content ah and spectator quark of flavour c 




: decuplet flavour wave functions 




: quark propagator of tlie spectator quark a 


D(ab) resp. Dg^J 


: scalar resp. axial vector diquark propagator 




: exchange kernel for quark flavour a, see eq.(^ for the 


(p, A = 1 ... 5) 


definition of S 



Of course, different to eqs. (|[]^), different masses corresponding to the flavour content 



of quark and diquark are used in the numerators and denominators of the propagators. 

For octet states each flavour wave function with scalar diquarks is to be expanded in 
Dirac space according to eq. (p!7| ) with only the scalar amplitudes 5*1 and 5*2 and each one 
with axialvector diquarks according to the same equation with the 6 vector amplitudes 
Ai . . . Aq. For decuplet states only flavour wave functions with axialvector diquarks are 
considered, they have to be expanded as indicated in eq. (p4D. 

The Bethe-Salpeter equation for S baryons now reads: 



/ 


{us)s \ 












\ 






usYs 


= -9l 












(59) 
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ssYu ) 
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9s " 



_9a 7>'I/,5 
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{us)s \ 






us^s 
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'ssYu J 



By interchanging s ^ u one obtains immediately the equation for S baryons. 
Broken SU{'i) couples the symmetric A and the flavour singlet. We introduce the 



flavour wave functions Fi = -^{(us)d — (ds] 
and the equation for the physical A reads: 



u 



F2 = {ud)s and A'^ = -^{[usfd - [ds] 



-u] 









^ SuD{us) 
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= -9l 
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-VI 


^9s « 










^ 9s " 9 


s 



I {us)s \ 
[usYs 

\\ssru) 

The S* baryons belonging to the decuplet are described by the equation: 



[ssYPu J V S, 



'uD^sd \V2K^'' )[[ssru 



(60) 



(61) 
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Eigenvalue vs. Bound State Mass 
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Figure 8: Eigenvalues vs. bound state mass M in the octet equation. Axialvector diquark 
propagator with and without diagonal approximation. 

As before the equation for the S* baryons may be obtained by interchanging s ^ u. 
Note that we neglected contributions from a state with scalar diquark (us). 

B Solutions with Proca Propagator for Axialvector 
Diquarks 

In order to study the influence of the full Proca propagator on vertex functions of spin-1/2 
baryons and the corresponding eigenvalues we are enforced to use a diquark size factor of 
the dipole type, writing instead of eq. (^) 



to regularise the equation which has no stable solution otherwise. 

In figure |^ we compare the eigenvalues as a function of M for both the Proca propagator 
and its diagonal approximation. Even in regions of moderate binding the eigenvalues do 
not differ by more than 2% which makes the approximation of the diagonal propagator 
in computing the masses a reliable one. While the vertex functions are essentially the 
same for both choices of the axialvector diquark propagator, the Ai components of the 
Bethe-Salpeter wave function differ by approximately a factor of 10. 




(62) 
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Figure 9: Eigenvalues vs. bound state mass M for direct and modified ladder approximation. 



C Direct vs. Modified Ladder Approximation 

In figure ^ we display for a representative parameter set the eigenvalues obtained in 
direct and modified ladder approximation. One clearly sees that for bound state masses 
approximately equal to the sum of the constituent masses the eigenvalues are almost 
identical, and even for strongly bound states the deviation is small. 



However, larger deviations occur in the vertex functions. Figure |10] shows zeroth and 
first Chebyshev moment of 5*1 for modified and direct ladder approximation. Whereas 
Si hardly differs for the two approaches, 5*^ receives a sign flip when switching between 
the two approaches. This occurs for almost all amplitudes in odd Chebyshev moments. 
When considering the Bethe-Salpeter equation for nucleons with 5'a=0, the different sign 
in S\ causes the electric form factor of the neutron to be changed. Isospin breaking 
effects are also different: The neutron-proton mass difference is in the direct approach 
bigger than the constituent quark mass difference — m„ whereas in the modifled ladder 
approximation it assumes values of approximately 0.7(mrf — iriu) [pOl . 
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Vertex Amplitude S, (p) 



Vertex Amplitude S, (p) 




modified ladder appr. 
■ direct ladder appr. 



2.0 3.0 

p/m^ 



4.0 5.0 6.0 



Figure 10: Zeroth (left) and first (right) Chebyshev moment of the vertex amplitude 5*1. 
Parameters are mq=mo+=mi+ , M=1.9mq, A=2mq, ga/gs=0-^ and d=10. 
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